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A dynamo model for axisymmetric and non-axisymmetric 
solar magnetic fields 
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ABSTRACT 

Increasing observations are becoming available about a relatively weak, but persistent, 
non-axisymmetric magnetic field co-existing with the dominant axisymmetric field on 
the Sun. It indicates that the non-axisymmetric magnetic field plays an important 
role in the origin of solar activity. A linear non-axisymmetric a 2 — Q, dynamo model 
is set up to discuss the characteristics of the axisymmetric (m = 0) and the first non- 
axisymmetric (m =1) modes and to provide further the theoretical bases to explain 
the 'active longitude', 'flip-flop' and other non-axisymmetric phenomena. The model 
consists of a updated solar internal differential rotation, a turbulent diffusivity varied 
with depth and an a-effect working at the tachocline in rotating spherical systems. 
The difference between the a 2 — and the a — fl models and the conditions to favor 
the non-axisymmetric modes with the solar-like parameters are also presented. 

Key words: magnetic fields - MHD - Sun: activity - Sun: magnetic fields. 



1 INTRODUCTION 



The distribution of magnetic field emerging on the solar 
surface carries clues to the mechanism of the field gen- 
eration. One striking feature of this distribution is clus- 
tering of active regions w hi ch is commonly called 'ac- 
tive longitude' dBai I fl 987; Bcn evolenskava et al. I 1 19991 : 
|Pe Toma, White fc Harvey I [2000). Signatures of possible 
longitudinal inhomogeneities have also been reported in 
the distribu tions of solar wind and interplanetary mag- 
netic field (|Neugebauer et all I2000T ) . Furthermore, 'flip- 
flop' phenomenon, i.e. the two persistent active longitudes 
separated by 180°, has also been identified on the Sun 
(|Berdvugina fc Usoskm~ll2003l) . These observations indicate 
the involvement of large-scale non-axisymmetric magnetic 
field in the formation and evolution of the dominant ax- 
isymmetric solar activities. Hence, it is valuable to set up 
the non-axisymmetric dynamo model to explain these non- 
axisymmetric solar magnetic fields. 

The pioneer works on the theoretical investigations of 
the non-axisymmetric activities can be classed mainly as 
two kinds. One is that the generation sources are non- 
axisymmetric, and the non-axisymmetric magnetic field is 
produced accordingly. For example, Bigazzi & Ruzmaikin 
(2004) and Moss et al. (2002) adopted the non-axisymmetric 
distribution of cv-effect. The other is based on the ax- 
isymmetric sources of generation but to excite the non- 
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axisymmetric field. The numerical results of Chan et al. 
(2004) supported this possibility. 

Earlier studies (|Stix I Il97ll : llvanova fc Ruzmaikin I 

1985) concerning the linear non-axisymmetric solar dynamo 
with decoupled axisymmetric and non-axisymmetric modes 
have been taken. But these earlier studies could not include 
the correct distribution of solar differential rotation, which 
was unknown at that time. Recently, there are some works 
on the non-linear non-axisymmetric dynamo models. Moss 
(1999) obtained stable solutions which possessed a small 
non-axisymmetric field component co-existing with a domi- 
nant axisymmetric part with the updated solar rotation pro- 
file. Bigazzi & Ruzmaikin (2004) studied the generation of 
non-axisymmetric fields and their coupling with the axisym- 
metric solar magnetic field. Bassom et al. (2005) used an 
asymptotic WKBJ method to investigate a linear a 2 — 
model with the aim to isolate the basic physical effects lead- 
ing to the preferable excitation of non-axisymmetric solar 
and stellar magnetic structure. However, are there possibil- 
ities to work out a linear non-axisymmetric solar dynamo 
with the updated generation sources to relate with the non- 
axisymmetric phenomena? What are the differences, such 
as configuration and cycle and so on, between the axisym- 
metric and non-axisymmetric modes? When will the non- 
axisymmetric mode be preferred? These are the main objec- 
tives of the paper. 

With the axisymmetric sources of generation, we de- 
velop a new high-precision non-axisymmetric code based 
on the spectral method and begin with the linear non- 
axisymmetric mean field dynamo equations. The axisym- 
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metric mode m = and the first non-axisymmetric mode 
m — 1 are discussed, respectively, in Sec. 5 and Sec. 6. We 
will show the difference between the a 2 — f2 and the a — fi 
models in Sec. 3. In Sec. 4, the condition to excite the dom- 
inant axisymmetric mode and the condition to favor the 
non-axisymmetric mode will be discussed. 



2 MATHEMATIC FORMULATIONS 
2.1 The basic equations 

The starting point of our model is the mean field dynamo 
equation, governing the evolution of the large-scale magnetic 
field B in response to the flow field U, the a-effect and the 
magnetic diffusivity rj: 

^ =V x[UxB + aB -?jV x B\. (1) 

Since the turbulent diffusivity is much larger than the molec- 
ular diffusivity, we ignore the molecular diffusivity in r\. For 
the flow field, only the (differential) rotation Q is considered 
for simplicity. Since the magnetic field is divergence-free, we 
expand B in term of two scalar functions h and g which rep- 
resent the poloidal and toroidal potentials , respectively, in 
the sp heric al polar c oordi nates (r, 8,<j>) as IChandrasekhar I 



trie sp heric al polar coordin, 1 
jl96lh and lMoffattl l|l978h . 

B = V x V x rh(r, 9, <j>, t) + V X rg(r, 8, <f>, t). 



(2) 



When a — a(r, 8), Q — Q(r, 8), r\ — rj(r), substituting equa- 
tion^) in equation (1), the governing equation reduces to: 



8 t L 2 h = R a V° + V V 2 L 2 h + RqV£ n + R a VZ r 

= R a V Q h + R n V& + R n Vf w + R a V° N 

t-,2 t 2 dt] d 2 1 dri 2 
+ L g+ —-L g+ -—L g, 

or or r or 



d t L 2 g 



where 



(3) 



(4) 



;( sm ^) 



i 8 o d 



V = --^-r -L 



r ■ V x (aB) 
r • V x (Ux B) 
r ■ V x [V x (aB)] 
Vx[Vx([/xB)] 



sin 8 89 y 89' sin 2 8 ' r dr dr 

Vqn and V a N are the terms which have the azimuthal com- 
ponent 8 /8(f> and Vh, V a ,VnN and V a N can be obtained 
from (see Appendix): 

Vl + V^N, (5) 

(6) 
(7) 
(8) 

The equations (3) and (4) have been cast in non-dimensional 
form by expressing all lengths in units of solar radius Rq and 
time in units of the magnetic diffusion time Rq/t/o. This has 
led to the appearance of two dimensionless numbers: 

Q-oRq 



V' 



+ V 9 



V£ + V£ N . 



Rat — 



Rn 



0) 



(10) 



where a a and r\ are reference values for the a-effect and the 
diffusivity in the convective zone (CZ), respectively. And 
f2 is the characteristic value of the differential rotation. 
The quantities R a and Rn are dynamo numbers measuring 
the relative importance of inductive versus diffusive effects. 
More discussion about -Rn will be given in Subsection 2.2. 
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Figure 1. Radial distributions of solar rotation at different lat- 
itudes. There are strong radial shear in the high latitude of the 
tachocline. At 35° latitude, the shear is very weak (dashed line). 



2.2 Internal rotation fi(r, 8) 

Based on the helioseismic inversion (|Schou et al. I Il998l : 
ICharbonneau et al. Ilf999l ). there are two strong radial shear 
regions inside of the Sun. One is in the tachocline and the 
other locates in sub-photospheric layer. For the sake of sim- 
plicity on computational solutions, We neglect the shear at 
the sub-surface and regard that the dynamo works in the 
tachocline. The following expression for the solar interior 
rotation is adopted. 

Q(r,9) =fi c + -[l + er/(2- 



d 



:)](n,(0)-n c ), 



(ii) 



where Q s (9) — Qbq + a-2 cos 2 9 + cm. cos 4 9 is the surface lati- 
tudinal rotation and 9 is co-latitude. The parametric values 
are set as r c = 0.7 Rq, d = 0.057?©, Q c /2n = 430.0 nHz, 
Q.eq/2-k = 455.8 nHz, a 2 /27r = -51.2 nHz, a 4 /27r = -84.0 
nHz. Fig. f shows the radial distribution of fl(r, 9) at differ- 
ent latitudes. It reveals that £1 depends weakly on depth in 
bulk of CZ. But in the tachocline, the rotation rate changes 
from almost uniform in the radiative interior to depth de- 
pendent in the CZ. Within the tachocline, rotation increases 
with distance from the core at low latitudes, while it de- 
creases at high latitudes. At intermediate latitudes (near 
35°, dashed line in Fig. 1) rotation is almost independent 
on the depth. 

Furthermore, we base our model on the rotating spher- 
ical systems with the rotation velocity Q c of the inner core. 
Thus the differential rotation in the rotating frame Q c is 

0,'(r,8) = ~[1 + er/(2^-^)](27r x 25.8) x 



(I. - 1.98 cos 2 8 - 3.26 cos 4 8) (nHz) 



(12) 



The differential rotation of the surface at the equator is (27rx 
25.8) nHz and we regard it as the characteristic value of the 
differential rotation f2 in Eq. (10). Hence, the value of Rn 



is adopted as 



8 x IP 10 mV 

Vo 



which is only decided by the 



reference value of the diffusivity r? . 



2.3 The diffusivity profile rj(r) 

We use the analytical expression of Dikpati & Charbonneau 
(1999) for the diffusivity profile as 



V(r)=r 1c +\\l + erf(2 r —^)], 



(13) 
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[h(r, 6, </>, t),g(r, 9, </>, t)} = [h(r, 9, 0), g(r, 9, 0)]e s 



(15) 
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Figure 2. Radial distributions of the a-effect (solid line) which 
mainly concentrate in the tachocline and the magnetic diffusivity 
r\ (dashed line). Values of y-axis are not scaled. 



whose distribution can be seen in Fig. 2 (dashed line). The 
diffusivity r\ in CZ, is dominated by the turbulence. In 
the stably stratified core, the diffusivity r) c is much lower 
because of the much less turbulence. In what follows we 
take ?7c/»7o=0.01. The transition from high to low diffusiv- 
ity occurs near the tachocline, which is coincident with the 
rotational shear layer. Here, r\ is far less definite and is 
widely known to fall in the range from 2 x 10 10 cm 2 s _1 to 
2 x 10 12 cmV 1 . 



2.4 The a-effect ce(r,9) 

The a-effect cannot yet be determined from observations. 
The dominated physical mechanisms responsible for it can 
be categorized as the following three types. (1) It works 
at the surface produced by the decaying of active regions 
jBabcock|[T96]] ; iLeighton I Il969l ) . (2) It is direc tly related 
to turbulent convective motions (Parker 1955). It exists 
throughout the whole CZ and changes sign near the bot - 
tom of the CZ jKrivodubskii Ill998l ; Kuzanvan et al. Il2003h . 
(3) It works at the tac hocline induced by t he hydrodynam- 
ical she ar instabilities (IDikpati et al. II2001T ) or MHD insta- 
bilities dThelen Il2000l ). It is possible that all of them work 
simultaneously inside of the Sun. Here, we only consider the 
a-effect concentrating in the tachocline with the following 
expression 



a(r,9) = Qo I[l + er/(2^i)] 

1 -[l-erf{2 r —^)] COS 9, (14) 

where n = O.675i? , r 2 = O.725i? . The solid line in Fig. 2 
shows the variation of a(r, 9) with r, which is mainly concen- 
trated in the tachocline. The common angular dependence 
cost? is adopted, the simplest guaranteeing antisymmetry 
across the equator. Moreover, we do not consider the a- 
quenching since only the linear solutions are sought. 



2.5 The numerical scheme 

Since the governing equations (3) and (4) are two coupled, 
linear homogeneous equations in h and g, with the given 
boundary conditions, we can look for their eigensolutions 
with the form 



where s is the eigenvalue and can be written as s = a + iu. 
Only the solution that neither grows nor decays (a ~ 0), 
i.e. the onset of dynamo actions is considered. The corre- 
sponding R a is the critical R a . The solution with the lowest 
dynamo number is the easiest to excite and is the most sta- 
ble one. In what immediately follows, we solve the dynamo 
equations numerically using spectral (Chebyshev-r) method 
(see Jiang & Wang 2006 for detail). 

Different azimuthal modes m are decoupled in linear 
theory. For given m, we expand h and g at the onset of the 
dynamo action in terms of Chebyshev polynomial T n (r) and 
surface harmonics p^e 1 ™ 1 '* in the meridional circular sector 
r € [0.6, 1.0], 9 e [0.0, 7r] as follows: 



n — l — m 
N L 

ff = EE <i T -( ar - b)Pr(cos9)e^, 

n — l—m 



(16) 
(17) 



where ar — b £ [—1, +1]. Here, TV and L are the truncations 
needed to get convergence. It varies with different dynamo 
number and different Q, tj, a profiles, cjj ; and c 3 n l are eigen- 
vectors. 

At two interface r = r = 1.0 and r = r; = 0.6, both 
magnetic field and the tangential electric field must be con- 
tinuous. The exterior r > 1.0 is a vacuum and eigensolu- 
tions are matched to a potential field. The radiative core is 
assumed to behave as a perfect conductor. We may obtain 
(see Schubert & Zhang 2001 for detail): 

at r = ri = 0.6, ^ ^ c^-lfPrtcosfl) = 0, (18) 

n 1 

^^c^[a(-l) n+ Vr < + (-l) n ]Pr(cos^ = 0, (19) 

n I 

at r = r o = 1.0, ^ ^c^PPtcostf) = 0, (20) 

n 1 

^^c^ ! [an 2 r + (l + /)]P i m (cost9) = 0. (21) 

n I 

As pointed out by Ivanova & Ruzmaikin (1985), the 
system of Eqs. (3) and (4) may be decomposed into two 
subsystems, i.e. odd or even parity with respect to the 
equatorial plane. We denote them by A and 5*. With the 
parameters adopted in our model, both odd and even 
parity solutions have nearly the same dynamo number 
for given mode m. Therefore, we cannot decide which 
kind of symmetric solution is excited easier. However, ac- 
cording to the observations and theore tical computations 
l|Moss 1120041 ; iFluri fc Berdvugirial 120041 '). the dipolar mode 
(AO) for axisymmetric field and the perpendicular dipo- 
lar mode (SI) for non-axisymmetric field are definitely 
identified on the Sun although there are some possibilities 
for some other non-axisymmetric modes exist on the Su n 
(|De Toma. White fc Harvevl l2000l : ISong fc Waiigl 120051 ). 
Hence we will only choose the two modes AO and 51 as 
the representatives of axisymmetric and non-axisymmetric 
modes, respectively and investigate AO and SI in detail be- 
low. We firstly discuss the difference between the a 2 — Q, 
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10 2 10 3 10 4 

Rn 

Figure 3. A critical R a -RQ plot for the a — Q (solid line) and 
the a 2 — Q (dashed line) models for the axisymmetric mode AO. 
When Rq ^3x 10 3 , the a — f2 model has larger critical R a than 
the a 2 — fl one. When Rq > 3 X 10 3 , the two models have the 
nearly same critical R a . For the a — Q. model, R Rq is about 
3840. 
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Figure 5. A critical R a -Ra plot for the a 2 — Cl model with the 
modes AO (solid line) and SI (dashed line). When Rq ^3x 10 3 , 
the mode 51 has smaller R a . It is contrary for Rq > 3 X 10 3 that 
the mode AO owns smaller R a . 




Figure 4. The ratio of the magnetic energy between the toroidal 
and poloidal components vs. Rq for the a — Q (solid line) and 
the a 2 — Q (dashed line) models for the mode AO. 



and a — Cl models and then give the condition for the Sun 
to excite the preferred non-axisymmetric mode. 



3 THE a 2 - Cl DYNAMO MODEL VERSUS THE 

a - Cl ONE 

For the mean field dynamo theory, poloidal field is created 
from toroidal field by the a-effect and toroidal field from 
poloidal field by two ways, i.e. the differential rotation (Cl- 
effect) and a-effect. The model including all the ingredients 
is called the a 2 — CI model. When Rq 3> R%, the a-effect 
as the toroidal s ource can be ignored. The a — Cl model is 
always adopted (|Zeldovich et al. II 19831 ') . Is the simple a — Cl 
model fine for the Sun? What quantitative conditions does 
it need to satisfy? What are the definite differences between 
the two models with solar-like parameters? 

We first discuss them based on the axisymmetric mode 
AO and enlarge the range of r/ from 8 x 10 10 cm 2 s _1 to 
8 x 10 12 cm 2 s . Thus we obtain Rq ranging from 10 2 to 
10 4 . For the a — CI model, the condition for the generation of 
undamped magnetic field is o nly determined by D — R a Rn 
|lvanova fc Ruzmaikin I Il985), Accordingly we obtain the 
straight (solid) line in Fig. 3 with logarithmic abscissa and 
D = R a RQ is about 3840. But for the a 2 — CI model, it is 



more complicated with the generation of toroidal field by 
the a-effect (the dashed line in Fig. 3). Comparing the two 
lines of Fig. 3, we can see that for the a 2 — Cl model, dynamo 
action is increased contrasting with the a — Cl model by the 
reduction of the critical R a when Rq is small (< 3 x 10 3 ). 
With the increasing of Rq , the difference for the correspond- 
ing critical R a between the two models decreases. When 
Rn = 3x 10 3 , the agreement between the two models reaches 
the level of 0.3%. 

Since the absolute scale for the strength of magnetic 
field is undermined by linear eigenvalues calculations, we 
define the ratio of the magnetic energy between the toroidal 
and poloidal components as (jCharbonneau fc Macgregorl 
l200ll ): 



e 



J B\ dV 



(22) 



where Bt 
d 2 h 



L 2 h . 



Adh 
r 06 



—dg/d0 and Bp 

eg for the axisymmetric model. Fig. 4 gives the en- 
araa 

ergy ratios between the toroidal and poloidal fields at the 
onset state for the two models with different Rq. When Rq 
is less than 3 x 10 3 , the a — Cl model produces smaller Et / Ep 
than the a 2 — Cl model. The larger Rq is, the less difference 
the two models have. The agreement between the two mod- 
els reaches 0.46% with Rq = 3 x 10 3 . When Rq > 3 x 10 3 , 
the energy ratios are closely in conformity with each other. 
Thus we can replace the a 2 — Cl model by the a — Cl model 
and the corresponding turbulent diffusivity r\ should be less 
than 2.67 x 10 11 cm 2 s"\ 

Let us see the non-axisymmetric mode SI simply now. 
Table 1 gives the difference for the two models with differ- 
ent Rq. When Rq > 4 x 10 3 {rj a < 2 x 10 11 cm 2 s _1 ), the 
difference between the two kinds of models is 0.73% and the 
two models can be regarded as the same. In a word, if the 
non-axisymmetric mode m = 1 is considered in the model, 
it is necessary for Rq to be larger than 4000, namely r\ is 
less than 2 x 10 11 cm 2 s _1 so that the a — Cl model is at the 
limit of a 2 — Cl model. It is fully reasonable to adopt the 
a — Cl model to replace the a 2 — Cl one. 
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Table 1. The critical R a for the a 2 — Q and the a — SI models based on the first non-axisymmetric mode SI. d = (a2 — al)/a2. al is 
the critical R a for the a — f2 model and al for the a 2 — Q, model. 



Rn 


500 


1000 


2000 


3000 


4000 


5000 


6000 


7000 


8000 


al 


4.92 


2.45 


1.545 


1.320 


1.239 


1.203 


1.190 


1.180 


1.179 


a2 


4.17 


2.33 


1.510 


1.306 


1.230 


1.198 


1.186 


1.177 


1.176 


d 


18.0 


5.15 


2.32 


1.07 


0.73 


0.42 


0.34 


0.25 


0.25 



4 AXISYMMETRIC VERSUS 

NON-AXISYMMETRIC MODE 

It is commonly regarded that the strong differential rota- 
tion works in favor of the axisymmetric mode. Without 
it, all single main-sequence stars with outer CZs have the 
non-a xisymmetric field configurations l)Riidiger fc Elstner I 
1 19941 ). The non- uniform rotation produces the observed 
dominant oscillatory dipolar field on the Sun. However, is 
the differential rotation expressed by Eq. (12) is enough to 
produce the dominant axisymmetric field? What condition 
does it need to satisfy to make the non-axisymmetric mode 
preferred? 

With the given rotation profile (12) of the Sun, we 
range Rn from 500 to 8000. Accordingly, r\ changes from 
1.0 x 10 11 cm 2 s _1 to 1.6 x 10 12 cm 2 s _1 . Fig. 5 displays the 
critical R a with different Rn for the a 2 — Q model. The ax- 
isymmetric mode AO (solid line) has lower critical R a when 
Rn > 3000 and will be preferred to excite. It is contrary 
for Rn < 3000 that the non-axisymmetric mode 51 will 
have lower critical R a and will be the preferred mode. The 
smaller Rn is, the easier 51 is excited. According to Sec. 3, 
it is not at the a — Q limit (Rn < 3000) when 5*1 is the 
preferred mode. In other words, it is impossible to favor the 
non-axisymmetric mode at the a — fl limit. To get the pre- 
ferred non-axisymmetric modes, the contribution of a-effect 
to generation of th e toroidal field cannot completely van- 
ish (|Radler Ill986al ). The farther it deviates from the limit, 
the more important roles the a-effect plays to produce the 
toroidal field and the easier the non-axisymmetric mode to 
excite. This is also consistent with the analytical results of 
Bassom et al. (2005) . Furthermore, they gave the reason that 
the wind-up of non-axisymmetric structures can be compen- 
sated by phase mixing inherent to the a 2 — Q dynamo. 

Moreover, rather than the differential rotation, Rn is 
the decisive parameter to decide which kind of mode is pre- 
ferred. We may also say that the turbulent diffusivity r\ is 
the key parameter since the differentia l rotation has been ba - 
sically determined by the observation (|jiang fc Wand[2007h . 
For the a 2 - D. model, when Rn < 3000, i.e. r\ > 2.67 x 
10 11 cm 2 s _1 the non-axisymmetric mode will be preferred. 

In the coming two sections, we will take the a 2 — Q, 
model with r) =l-6 x 10 11 cm 2 s~' and _Rn=5000 so that the 
axisymmetric mode will be preferred. This is the real picture 
of the Sun. 
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Figure 6. Butterfly diagram of the toroidal field for the mode 
AO at the depth r = 0.7-Rq (left) and B r at the surface r = Rq 
(right). Solid (dashed) contours correspond to positive (negative) 
magnetic field. The diffusion time Rq/t) is taken as the time 
unit. 
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Figure 7. Evolution of the toroidal magnetic field in a meridional 
plane = separated by l/6th of the dynamo period T for 
the axisymmetric mode AO. Solid (dashed) contours correspond 
to positive (negative) toroidal field. The dot-dashed lines locate 
at 0.7Rq. The fields concentrate in the tachocline and the high 
latitudes where the radial shear is strong. 



5 THE AXISYMMETRIC MODE 

Left part of Table 2 is the truncation levels in the cal- 
culation of axisymmetric mode AO and the corresponding 
critical R a and frequency uj. N is the radial harmonics 



expanded in terms of Chebyshev function and L is the 
harmonics in Legendre function (see Eqs.(16)(17)). When 
N — 38 and L = 50, it gets convergence. The critical R a 
is 0.765 and frequency u> is ±410.37. With the dimension- 
less time t — Rq/t]o, we obtain the period is about 15 
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year, which is a bit shorter than the 22-years solar mag- 
netic cycle. In fact, the meridional circulation omitted in 
the p aper plays an important role in determining the true 
cycle l|Dikpati fc Charbonneau~lll999l ). The symbol '±' cor- 
responds to the dynamo wave propagating equatorward or 
poleward. 

Fig. 6 (left) displays the toroidal field at r = 0.7Rq with 
lu — +410.37. The magnetic field concentrates in the region 
from latitude 45° to 70° with strong radial shear dQ./dr < 
there. It is well-known that, without meridional circulation, 
the propagation direction of the dynamo wave is decided by 
the Parker- Yoshimura sign rule and it should be equator- 
ward since the product of the a-e ffect and radial gradient of 
differ ential rotation is negative (jParker 1 1 19551 : lYoshimura I 
1975). The solution with lj = —410.37 corresponds to the 
dynamo wave propagating poleward. Although some high- 
frequency dynam o waves were identified to prop agate pole- 
ward on the Sun (jMakarov II 19891 : lObridko et al. Il2006h . the 
poleward solution obtained in our method cannot be used 
to explain these observations. It is beyond the limitation 
of our method. Hence, the poleward solution is meaningless 
and should be neglected. Fig. 6 (right) is a time-latitude 
diagram of the radial field at the surface. The phase shift 
between the two comp onents is near i t/2 which is consistent 
with the observations l|Sheelev Ill99lh . 

Fig. 7 is the evolution of the toroidal field in a merid- 
ional plane <f> = at an interval of l/6th of solar cycle period. 
The dot-dashed line is 0.7Rq. The magnetic field is among 
the region where the radial shear of differential rotation is 
strong, i.e. the high latitude of the tachocline. This follows 
the general rule that differential rotation tends to destroy 
any deviation from axisymmetry and toroida l field favors to 
be produced in the strong ra dial shear region l|Moffatt II 19781 : 
iBieazzi fc Ruzmaikin Il2004 ). 

It seems that the location of the toroidal field pro- 
duced in the model is higher than that of the obser- 
vation. Furthermore, when the toroidal flux rope rises 
through the CZ to emerge, it will have the poleward 
deflection further due to the effect of Coriolis force 
l|Caligari. Moreno-Insertis. Schiissler I |l995). But in the 
model, we omit an important ingredient, i.e. the meridional 
circulation. If a meridional circulation is considered, the 
strong field produced within the tachocline at high latitude 
will be carried to the low latitudes. The toroidal flux entering 
the CZ will become buoyantly unstabl e and emerge to form 
the a ctive regions at the low latitudes (|Nandv fc choudhuri I 
l2002h . 



6 THE NON-AXISYMMETRIC MODE 

Right part of Table 2 is the truncation levels of the mode 
SI. When N = 28 and L = 70, it gets convergence, which 
is slower than the calculation for the axisymmetric one. 
The critical R a is 1.203, much larger than that for the ax- 
isymmetric mode. Therefore with taking Rn = 5000 and 
rj = 1.6 x 10 11 cm 2 s _1 , the axisymmetric mode will be 
preferred. Fig. 8 is the near surface distribution of the ra- 
dial magnetic field for the mode SI at the fixed time. The 
toroidal magnetic field of the mode m = 1 superimposed on 
the axisymmetric toroidal field produces a localized maxi- 
mum ('hump'). The non-axisymmetric enhancement of the 
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Figure 8. Near surface distribution of the radial component for 
the mode 51 at a fixed time. Solid (dashed) contours correspond 
to positive (negative) magnetic field. 



underlying magnetic field causes the cluster ing of sunspots 
to form 'active longitudes' l|Ruzmaikin Il200lh and 'Flip-flop' 
which behaves as a special phenomenon of 'active longi- 
tudes'. 

The frequency ui is -1.022. Hence its period is nearly 
400 times longer than that of the mode AO. Thus the 
non-axisymmetric mode SI appears to be rather steady or 
weakly oscillating c omparing to the axisymmetric mode AO 
(|Berdvugina I I2004T ) . The time variations of the mode AO 
are periodical. By changing the sign of the mode AO, the 
predominant longitude jumps by about 180°, which is just 
the flip-flop phenomenon (see the details of Sec. 3 of Fluri 
& Berdyugina (2004)). But only based on these two modes, 
the full flip-flop cycle has the same length as the AO cycle 
rather than the value which is a 3-4 times short er (about 3.7 
years for the Sun) than the main activity cycle l|Berdvugina I 
120041 ). There should have more complicated field configura- 
tion working on the Sun. 

Both the non-axisymmetric and axisymmetric mag- 
netic fields are generated by the same axisymmetric 
sources. They evolve independently on each other in lin- 
ear theory. In fact, some nonlinearities, such a s non- 
axisymmetr ic a-effect (Bigazzi fc Ruzmaikin~l l20o3 ). MHD 
instability dDikpati fc Gilman 20051 ). magne tic buoyancy 
JChatteriee et al. |2004 ) and the a- quenching (|Zhang et al. I 
120031 : iMoss 1 120051 ) induce the differe nt modes cou pled to- 
gether and produce the flip-flop cycle (|Moss 1120041 ). 

Left of Fig. 9 shows the contours of toroidal magnetic 
field in a meridional plane. All the fields concentrate in the 
lower part of the tachocline where the diffusivity is less than 
that of other regions (see Fig. 2). Right of Fig. 9 displays 
the butterfly diagram for the mode SI at the depth r = 
O.75i?0. The field is mainly concentrated around the 35° 
latitude, where the radial shear is weak (see Fig. 1). This 
is consistent with the work of Bigazzi & Ruzmaikin (2004) 
that the non-axisymmetric field survives only in the weak 
differential rotation and low diffusivity region. 



7 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the properties of the 
axisymmetric and the non-axisymmetric modes with a lin- 
ear a 2 — model in a rotating frame trying to understand 
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Table 2. Truncation levels and the corresponding critical R a and frequency u> for the axisymmetric mode AO (left) and the first 
non-axisymmetric mode SI (right). N and L are the harmonics in Chebyshev function and Legendre function respectively. 
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Ra 




34 


46 


0.762 


±408.97 


36 


46 


0.767 


±412.06 


38 


46 


0.764 


±409.90 


38 


48 


0.765 


±410.39 


38 


50 


0.765 


±410.37 




Figure 9. Contours of the toroidal field in a meridional plane 
(left) and Butterfly diagram of the toroidal field at the depth 
r = O.75.R0 (right) for the mode SI. The field concentrates in 
the low diffusivity (below the dot-dashed line 0.7Rq in the left 
diagram) and weak differential rotation (about the latitude 35° 
in the right diagram) region. Solid (dashed) contours correspond 
to positive (negative) magnetic field. 



the active longitudes, flip-flops and other non-axisymmetric 
phenomena. The model consists of updated differential rota- 
tion, turbulent diffusivity varied with depth and the a-effect 
working in the tachocline. The definite differences between 
the a — Q and a 2 — fl models and the conditions to excite the 
non-axisymmetric modes with the solar-like parameters are 
presented. The Chebyshev-r method is used to numerically 
solve the problem with a high precision. We have tested this 
code with the analytical solutions of a simple a 2 model and 
compared with the results of Stix (1976). They matched well 
with each other. 

Although the conditions to favor the non-axisy mmetric 
modes h ave been investigated by some works ijRadler I 
Il986al lbl; H assom. Kuzanvan. Sokoloff. Soward I 120051 ). we 
are the first to apply the updated solar parameters with 
numerical method to obtain these. We also point out the 
role of the turbulent diffusivity r\ during the process. It is 
poorly known and cannot be directly deduced from obser- 
vations. When it is lower than 2.0 x 10 11 cm 2 s _1 and Rn 
is larger than 4000 accordingly, the a — model is at the 
limit of the a 2 —Q, model. Furthermore, based on the a 2 — Q, 
model, when R a is lower than 3000, the non-axisymmetric 
mode m — 1 will be more and more easily excited than 
the axisymmetric mode m = when increasing the value of 
diffusivity r) and decreasing the corresponding Rq . 

The non-axisymmetric mode SI has much longer period 
than the axisymmetric one AO. The co-existence of oscillat- 



JV 
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Ra 




24 


66 


1.201 


-0.794 


26 


66 


1.198 


-0.919 


28 


66 


1.198 


-0.967 


28 


68 


1.203 


-0.946 


28 


70 


1.203 


-1.022 



ing mode AO and the nearly steady SI results in the flip- 
flop phenomenon. Because the differential rotation affects 
axisymmetric and non-axisymmetric magnetic fields in the 
different ways, the two kinds of fields prefer to occur in differ- 
ent regions. The axisymmetric magnetic field is mainly con- 
centrated near the high latitude (about 55°) around 0.7 Rq, 
where the radial shear of differential rotation is strong. How- 
ever, the non-axisymmetric field occurs near the intermedi- 
ate 35° latitude in the bottom of the tachocline, where the 
differential rotation is weak and the diffusivity is low. 

Usoskin et al. (2005) gave that the ratio between the 
non-axisymmetric strength and the axisymmetric one is 
roughly 1:10 by analyzing sunspot group data for the past 
120 years. Based on the non-linear models, Moss (1999, 
2004) presented the energy ratios between the two kinds of 
modes although the values cannot match well with the ob- 
servation. Since our model is linear and different modes are 
decoupled, it is the limitation to provide the energy ratios 
between the two modes. 

In this work, we regard the strong radial shear only 
exists in the tachocline and omit the sub-surface shear 
and other details of the distribution of differential rotation. 
Brandenburg (2005) argued for the alternative ideas con- 
cerning dynamo operating in the bulk of CZ, or perhaps 
even in the sub-surface shear layer. Moreover, we also only 
take the a-effect working in the tachocline and ignore the 
other two possible mechanisms. The two generation sources 
are still the hot topics on debate. Since we do not aim to 
give the detailed description of the Sun and just put em- 
phasis on the basic characters of the axisymmetric and non- 
axisymmetric modes, it is feasible for us to take the two 
simple generation sources and set up the thin-layer dynamo 
model. Of course, more rich and realistic models can open 
new option for the understanding of solar magnetic field. 

Furthermore, since we tried to expatiate on our objec- 
tives with the simple generation sources, the meridional cir- 
culation is not considered in the wor k. It plays an impor- 
tant role in the axisymmetric m ode (|Nandv fc choudhuri I 
120021 : iGuerrero fc Mufiozl |2004| ) . It carries the strong ax- 
isymmetric toroidal field produced at the high latitudes to 
the low ones and produces the active regions there with the 
magnetic buoyancy. But it has no much influence on the 
non-axisymmetric field according to the work of Bigazzi & 
Ruzmaikin (2004). 

In the forthcoming studies we will include the nonlin- 
earities and the meridional circulation to investigate the in- 
fluence on the coupling of the different modes and the role 
of meridional circulation in the non-axisymmetric dynamo. 



8 J. Jiang and J. X. Wang 



ACKNOWLEDGMENTS 

We thank the anonymous referee for valuable comments 
which have helped in improving the paper. We extend our 
thanks to A. R. Choudhuri and K. M. Kuzanyan for review- 
ing the manuscript and giving many helpful comments. JJ 
would like to thank X. H. Liao and K. K. Zhang for kind 
supervision on the work. This work has been supported by 
National Natural Science Foundation of China (10573025, 
10603008) and by the National Basic Research Program 
(G2006CB806303). 



REFERENCES 

Babcock H. W., 1961, ApJ, 133, 572 
Bai T., 1987, ApJ, 314, 795 

Bassom A. P., Kuzanyan K. M., Sokoloff D., Soward A. M., 
2005, Geophys. Astrophys. Fluid Dynam., 99, 309 

Benevolenskaya E. E., Hoeksema J. T., Kosovichev A. G, 
et al., 1999, ApJ, 517, L163 

Berdyugina S. V., Usoskin I. G, 2003, A&A, 405, 1121 

Berdyugina S. V., Moss D., Sokoloff D., Usoskin I. G, 2006, 
A&A, 445, 703 

Berdyugina S. V., 2004, Sol. Phy., 224, 123 

Bigazzi A., Ruzmaikin A., 2004, ApJ, 604, 944 

Brandenburg A., 2005, ApJ, 625, 539 

Caligari P., Moreno-Insertis F., Schiissler M., 1995, ApJ, 
441, 886 

Chan K. H., Liao X., Zhang, K., Jones C. A., 2004, A&A, 
423, L37 

Chandrasekhar S., 1961, Hydrodynamic and hydromag- 

netic stability. Oxford: Clarendon Press 
Charbonneau P., Christensen-Dalsgaard J., Henning R., 

Larsen R. M., Schou J., Thompson M. J., Tomczyk S., 

1999, ApJ, 527, 445 
Charbonneau P., Macgregor K. B., 2001, ApJ, 559, 1094 
Chatterjee P., Nandy D., Choudhuri A. R., 2004, A&A, 

427, 1019 

Choudhuri A. R., 2003, Sol. Phy., 215, 31 
De Toma G, White O. R., Harvey K. L., 2000, ApJ, 529, 
1101 

Dikpati M., Charbonneau P., 1999, ApJ, 518, 508 

Dikpati M., Gilman P. A., 2001, ApJ, 551, 536 

Dikpati M., Corbard T., Thompson M. J., Gilman P. A., 

2002, ApJ, 575, L41 

Dikpati M., Gilman P. A., 2005, ApJ, 635, L193 
Fluri D. M., Berdyugina S. V., 2004, Sol. Phy., 224, 153 
Guerrero G. A., Munoz J. D., 2004, MNRAS, 360, 317 
IvanovaT. S., Ruzmaikin A. A., 1985, Astron. Nachr., 306, 
177 

Jiang J., Wang J. X., 2006, Chin. J. Astron. Astrophys., 6, 
227 

Jiang J., Wang J. X., 2007, J. Adv. Space Res., 

doi:10.1016/j.asr.2007.01.033 
Krivodubskii, V. N., 1998, Astronomy Reports, 42, 122 
Kuzanyan K. M., Lamburt V. G, Zhang H. Q., Bao S. D., 

2003, Chin. J. Astron. Astrophys., 3, 257 
Leighton R. B., 1969, ApJ, 156, LI 

Obridko V. N., Sokoloff D. D., Kuzanyan K. M., Shelting 
B. D, Zakharov V. G, 2006, MNRAS, 365, 827 
Makarov V. I., 1989, Sol. Phy., 123, 367 



Moffatt H. K. 1978, Manetic Field Generation in Electrical 

Conducting Fluids. Cambridge Univ. Press, Cambridge 
Moss D., 1999, MNRAS, 306, 300 
Moss D., 2004, MNRAS, 352, L17 
Moss D., 2005, A&A, 432, 249 

Moss D., Piskunov N., Sokoloff D, 2002, A&A, 396, 885 
Nandy D., Choudhuri A. R., 2002, Science, 296, 1671 
Neugebauer M., Smith E. J., Feynman J., Feynman J., 

Vaughan A. H., 2000, J. Geophys. Res., 105, 2315 
Parker E. N., 1955, ApJ, 122, 293 

Radler K. H., 1986, ESA SP-251, Plasma Astrophysics, 

ESA Publications Division, Noordwijk, 569 
Radler K. H., 1986, Astron. Nachr., 307, 89 
Riidiger G, Elstner D, 1994, A&A, 281, 46 
Ruzmaikin A., 2001, Space Science Reviews, 95, 43 
Schou J., Antia H. M., Basu S. et al., 1998, ApJ, 505, 390 
Schubert G, Zhang K. K., 2001, ApJ, 557, 930 
Sheeley N. R. Jr, 1991, ApJ, 374, 268 
Song W. B., Wang J. X., 2005, ApJ, 624, L137 
Stix M., 1971, A&A, 13, 203 
Stix M., 1976, A&A, 47, 243 
Thelen J.-C, 2000, MNRAS, 315, 165 
Usoskin I. G, Berdyugina S. V., Poutanen J., 2005, A&A, 
441, 347 

Yoshimura H., 1975, ApJ, 201, 740 

Zeldovich Ya., Ruzamikin A., Sokoloff D., 1983, Magnetic 
fields in Astrophysics, Gordon and Breach, New York. 

Zhang K., Chan K. H., Zou J., Liao X., Schubert G, 2003, 
ApJ, 509, 663 



APPENDIX A: THE FULL FORM OF THE 
GOVERNING EQUATIONS 

The full forms of the governing equations about the toroidal 
field g and poloidal field h in (3) and (4) are as follows: 



dL 2 h 
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0L 2 g 
dt 



= R a [aL g 



dadg T7 2 r 2, 



Rn[n^] + R a 
R a [-aV 2 L 2 h + 



1 daAdh d 2 h 



sin 6 88 r d<f> drd(f> 



(Al) 



+ 



drd6 drd8 



h 



+ 



da dL 2 h _ 
dr dr 
1 da dL 2 h 

'V 2 ~m 39 

1 i 



Ida 2 
-— L h- 

r or 

i d 2 h 



9 , . „ da 2 , 
— (sine — )L h 



r drd6 



r sm do dO dr do 

i d dn d dh 

[sm6 frd6 {sme d6 )] 



suit 
1 dQ & 



dh 



+ 



sin 9 dr d(j) 2 d9 r sin 6 d(j) 2 dr 

„2 r 2 . dr/ d 2 . 1 dr) 2 
r)V L g+TTjrL g+- — L g 
dr dr r dr 



d an 



+ r f 1 dff , d a Ida dg 
a[ sm6d<f> ( drd6 rd6 dr 1 



(A2) 



Non-axisymmetric solar dynamo 

where a, Q and rj are the expressions after non- 
dimensionalization. 

This paper has been typeset from a Tj^X/ I^Tf^X file prepared 
by the author. 



